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Section  1 

INTRODUCTION  AND  BACKGROUND 


Vertical  shear  is  the  vertical  gradient  of  hori¬ 
zontal  velocity,  and  is  important  from  several  points  of 
view.  Details  of  the  Navy  interest  in  shear  were  summarized 
in  an  SAI  report  to  NORDA  (Lambert  and  Patterson;  1980), 
which  reviewed  the  status  of  measurements  of  vertical  shear 
in  the  upper  ocean  as  of  early  1980.  The  results  of  that 
report  led  us  to  focus  our  attention  on  describing  and 
modeling  shear  over  the  scale  range  between  1  meter  and  a  few 
tens  of  meters. 

Navy  applications  include  the  potential  effects  of 
shear  on  subsurface  wakes.  Shear  can  kinematically  distort  a 
passive  wake  contaminant,  and  can  interact  with  active  wake 
signature  fields,  such  as  turbulence  and  internal  waves.  The 
oceanographic  community  is  generally  interested  in  fine-scale 
shear  distribution  because  of  the  (still  controversial) 
importance  of  shear  instability  in  turbulence  and  mixing. 
The  first  approaches  evaluated  the  contribution  of  internal 
wave  motions  to  the  shear  field  (Grabowski;  1980),  and  these 
focused  first  on  vertical  profiles  and  distributions  in  terms 
of  vertical  wavenumbers.  These  studies  led  to  the  confir¬ 
mation  of  earlier  studies,  that  the  variability  in  the  shear 
field  is  concentrated  in  the  high-wavenumber  end  of  the 
spectrum  and  in  the  low  frequency  (near-inertial)  portion  of 
the  internal  wave  band. 

A  number  of  other  studies  have  focused  on  statis¬ 
tical  descriptions  of  the  probability  distribution  of  inter¬ 
nal  waves,  concentrating  on  velocity  or  temperature  fluc¬ 
tuations,  but  not  specifically  on  shear.  Briscoe  (1977) 
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compiled  the  results  of  eight  such  studies.  He  noted  that 
the  most  impressive  result  of  his  compilation  is  that  con¬ 
tradictory  results  seem  to  emerge  from  similar  experiments. 
Comparisons  among  evaluations  of  the  Gaussianity  of  tem¬ 
perature  fluctuations  and  currents  yield  no  clear  consensus. 
Briscoe  suggested  that  this  may  indicate  intermittency  in  the 
structure  of  the  probability  distribution.  He  examined  the 
time  variability  of  horizontal  current,  temperature,  and 
vertical  displacement  data  from  the  Internal  Wave 
Experiment  (IWEX).  In  the  case  of  horizontal  currents, 
he  found  that  interludes  of  non-Gaussian  statistics  are 
related  to  kurtosis  values  of  less  than  3. 

To  the  authors'  knowledge,  only  two  studies  of 
shear  statistics  have  been  previously  performed,  both  using 
data  from  free  fall  vertical  profilers.  Again,  the  results 
were  contradictory.  Simpson's  (1975)  observations  of  shear 
probability  density  functions  over  a  length  scale  of  about 
35  cm  were  highly  intermittent  and  non-Gaussian.  Hebenstreit 
and  Grabowski  (1981)  found  that  shear  probability  density 
functions  over  2,  4,  8,  and  16  m  were  almost  all  nearly 

Gaussian . 

Vertical  profiles  of  shear  show  a  high  degree  of 
correlation  between  shear  squared  (S^)  and  Vaisala  fre¬ 
quency  squared  (N^).  Patterson  et  al .  (1981)  documented 

this  correlation  using  profiles  having  vertical  resolution 
between  2  m  and  8  m.  They  also  found  that  the  correlation 
improves  when  and  are  low-pass  filtered  in  the 

vertical.  The  implication  is  that  and  are  better 

correlated  over  length  scales  in  the  range  between  5  m  and 
20  m  than  over  a  2m  length  scale. 
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The  majority  of  vertical  shear  observations  have 
been  made  using  vertical  profiling  instruments.  Gargett  et 
al.  (1981)  describe  a  composite  vertical  wavenumber  spectrum 
of  upper  ocean  shear  from  three  different  velocity  profilers. 
Each  of  the  instruments  is  effective  over  a  limited  range  of 
vertical  length  scales,  and  by  carefully  comparing  their 
spectra,  Gargett  et  al.  form  a  composite  shear  spectrum  over 
vertical  wavelengths  ranging  from  a  few  centimeters  up  to 
more  than  100  m.  In  the  low  wavenumber  portion  0  <  k  <_  .1 
cpm,  they  find  that  the  shear  spectral  density  is  approx¬ 
imately  independent  of  wavenumber,  and,  as  a  consequence, 
scales  according  to  N2 .  This  result  is  consistent  with  the 
previously-noted  point-wise  correlation  between  shear  and 
s tratif ication.  The  results  are  most  reliable  in  the  deeper 
thermocline  where  N  does  not  vary  rapidly  with  depth. 

A  similar  correlation  is  observed  in  current  meter 
time  series--a  quite  different  type  of  data  set.  Eriksen 
(1978)  measured  velocity  differences  over  a  fixed  separation 
of  7  m  in  the  main  thermocline  near  Bermuda,  and  found  that 
the  shear  was  constrained  such  that  the  Richardson  number 
Ri  =  N2/S2  was  almost  always  greater  than  about  1/4. 
Davis  et  al.  (1981)  observed  a  very  similar  Richardson  number 
constraint  acting  on  15-min  averages  of  shear  calculated 
over  3  m  vertical  spacing  during  MILE  (Mixed  Layer 
Experiment)  . 

In  the  seasonal  thermocline,  N2  varies  rapidly 
with  depth.  Thus  shear  statistics  derived  from  vertical 
profiler  data  in  the  seasonal  thermocline  are  nonstationary 
in  depth.  Therefore,  vertical  wavenumber  spectra  of  shear 
are  probably  not  statistically  reliable,  especially  in  the 
low  wavenumber  band.  But  it  is  possible,  using  pairs  of 
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current  meters  separated  by  a  range  of  vertical  separations 
Az ,  to  extract  useful  information  about  the  vertical  scales 
of  shear  in  the  seasonal  thermocline.  One  of  the  chief 
objectives  of  this  report  is  to  evaluate  a  statistical  model 
which  is  capable  ot  relating  the  statistics  of  shear 
calculated  from  vertical  profilers  to  the  statistics  of  shear 
calculated  from  current  meters.  The  model  is  derived  from  an 
idealized  version  of  the  Gargett  et  al.  spectrum;,  and  the 
data  set  used  is  the  set  of  current  meter  records  from  MILE. 
The  success  of  this  model  allows  us  to  extend  the  range  of 
validity  of  the  composite  spectrum  of  Gargett  et  al.  (1981) 
into  the  seasonal  thermocline. 

This  report  is  organized  as  follows.  In  Section  2 
we  discuss  the  MILE  data  set,  including  the  current  meter 
response,  the  contamination  of  current  measurements  due  to 
mooring  motions,  and  the  low-pass  filter  used  to  remove 
unwanted  noise.  We  also  show  how  N2  was  computed  from  time 
series  of  AT/Az  combined  with  T-S  relations  calculated  from 
averaged  CTD  profiles.  In  Section  3  we  discuss  the  temporal 
variability  of  shear  distribution,  showing  that  shear  is 
dominated  by  inertial  oscillations. 

In  Section  4  we  describe  the  development  of  the 
statistical  model  of  mean-square  shear  based  on  an  idealized 
version  of  the  composite  vertical  shear  spectrum.  The  model 
predicts  mean  square  shear  as  a  function  of  vertical  sep¬ 
aration  Az  and  mean  square  Vaisala  frequency,  and  is  compared 
with  shear  observations  from  current  meter  time  series.  In 
Section  5  we  use  this  model  to  show  that  shear  distribution 
recorded  over  many  inertial  periods  is  reasonably  Gaussian. 

In  Section  6  we  consider  the  correlation  between 
N2  and  S2  over  very  short  time  scales  (15-minute 
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averages).  The  correlation  is  shown  to  be  significant,  and 
we  also  show  from  the  distribution  functions  of  s2/\2f 
that  knowing  the  value  of  allows  one  to  predict 

within  a  factor  of  two  at  least  46%  of  the  time,  and  within  a 
factor  of  four  at  least  75%  of  the  time. 

In  Section  7  we  sunmarize  the  conclusions  and 
recommend  possible  directions  for  future  work. 


Section  2 
DATA  PROCESSING 


2.1  CURRENT  METERS  AND  FILTERING 

The  current  meter  data  from  the  MILE-1  mooring 
provide  good  vertical  resolution  of  currents  over  a  19  day 
period  (19  August  to  6  September,  1977).  The  mooring  was 
deployed  in  the  vicinity  of  Ocean  Weather  Station  P,  at 
49‘37'N,  145'6'W,  where  the  ocean  depth  is  4360  m.  There 
were  19  vector  averag ing  current  meters  (VACM's)  deployed  at 
mean  depths  of  5,  8,  11,  14,  20,  23,  26,  29,  32,  35,  38,  41, 
44,  47,  50,  70,  92,  125,  and  175  m.  A  vector  measuring  cur¬ 
rent  meter  (VMCM)  was  located  at  94  m. 

The  VACM  records  current  and  temperature  at  112.5 
sec  intervals  (32  cph) .  The  VACM  measures  the  polar  coor¬ 
dinates  of  the  velocity  field  using  a  Savonius  rotor,  a  vane, 
and  a  magnetic  compass.  In  a  sense,  the  VACM  is  analogous  to 
an  odometer,  rather  than  a  speedometer.  Every  1/8  rotor  turn 
(corresponding  to  about  4.7  cm  of  wafer  displacement),  the 
sine  and  cosine  of  the  angle  of  the  vane  relative  to  magnetic 
north  are  accumulated  in  east-west  and  north-south  registers. 
The  registers  are  accumulated,  and  are  transferred  digitally 
to  a  magnetic  tape  (and  reset)  each  recording  interval. 
Later  in  the  laboratory,  the  records  are  processed  and  con¬ 
verted  into  physically  meaningful  measurements  (McCullough, 
1975)  . 

Halpern  et  al.  (1981)  performed  intercomparison 
tests  between  various  types  of  current  meters  and  moorings 
deployed  during  MILE  and  JASIN  (Joint  Air-Sea  Interaction 


Project).  They  compared  spectra  between  VACM's  and  VMCM's, 
and  among  surface  toroid,  spar,  and  torpedo  buoys,  and  sub¬ 
surface  buoy  moorings.  From  each  comparison  they  determined 
a  maximum  usable  frequency,  defined  as  the  frequency  beyond 
which  the  spectra  diverge.  The  maximum  usable  frequencies  of 
the  VACM's  deployed  in  the  upper  20  m  ranged  from  4  to  8  cph. 
In  deeper  water  the  maximum  usable  frequencies  are  smaller 
due  to  smaller  velocity  fluctuation  magnitudes.  During 
JASIN,  at  depths  between  79  m  and  92  m,  maximum  usable 
frequencies  ranged  from  0.2  to  8  cph.  During  MILE  only  a 
single  estimate  of  the  maximum  usable  frequency,  2  cph,  was 
available  in  this  depth  range. 

In  order  to  eliminate  the  high  frequency  oscil¬ 
lations  which  result  from  mooring  motions  and  imperfect  cur¬ 
rent  meter  response,  we  have  applied  a  low  pass  filter  to  the 
current  meter  records.  We  used  a  special  FORTRAN  program 
(Rabiner  and  Gold,  1975)  written  by  Jim  McClellan  at  Rice 
University  to  design  an  optimal  finite  impulse  response 
filter.  We  designed  the  128-point  symmetric  low  pass  filter 
to  have  a  cutoff  frequency  of  2  cph  and  a  transition  band 
from  2  cph  to  2.5  cph.  The  original  records  were  sampled  at 
32  hr*^-,  so  that  1  hr  of  data  was  lost  from  the  beginning 
and  the  end  of  the  record.  We  subsampled  the  filtered  record 
at  8  hr~l  (Nyquist  frequency  «  4  cph).  This  sampling 

rate  was  deliberately  chosen  to  be  a  little  higher  than 
necessary,  in  order  to  minimize  possible  aliasing. 

Since  we  are  interested  in  the  probability  density 
distributions  of  shear,  a  question  arises  concerning  the 
effects  of  instrument  nonlinearity  and  calibration  uncer¬ 
tainty  on  the  observed  distributions.  Briscoe  (1977)  con¬ 
sidered  a  simple  statistical  model  of  these  effects  on 
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VA CM's.  From  the  calibration  carves  of  McCullough  (1975),  he 
estimated  that  these  effects  will  adversely  affect  skewness 
and  kurtosis  observations  by  only  a  few  percent.  He  con¬ 
cluded  that  difficulties  with  VACM  calibration  uncertainty 
are  very  unlikely  to  be  the  source  of  observed  large  devi¬ 
ations  from  Gaussianity. 

2.2  DENSITY  STRATIFICATION  FROM  TEMPERATURE  GRADIENTS 

Figure  2.1  shows  temperature  and  salinity  profiles 
from  two  CTD  casts  on  21  and  24  August  during  MILE.  During 
this  period,  the  surface  mixed  layer  deepened,  but  the  main 
thermocline  shallowed  by  about  the  same  amount,  and  the  water 
column  above  100  m  has  cooled  and  become  more  saline. 

Figure  2.2  shows  the  averaged  temperature- 
salinity  relations  from  CTD  casts.  Above  75  m  (below 
ot  =  26  gm/cm3)  ,  about  80%  of  the  density  gradient  is  due 
to  the  temperature  gradient,  and  20%  is  due  to  the  salinity 
gradient.  Between  100  and  150  m,  virtually  all  of  the 
density  gradient  is  due  to  salinity. 

Temperature  sensors  were  potted  within  the  VACM 
housings.  Because  of  the  thermal  inertia  of  the  housings, 
the  temperature  response  is  attenuated  at  frequencies  above 
1  cph  (Levine,  1980).  In  order  to  produce  temperature  rec¬ 
ords  with  the  same  sampling  intervals  as  those  for  current, 
we  applied  the  same  filtering  process  described  in  Section 
2.1.  There  were  no  conductivity  sensors  attached  to  the 
current  meters.  Estimates  of  time-dependent  density  strati¬ 
fication  are  thus  derived  from  the  temperature  sensor  records 
alone,  using  mean  salinity  profiles,  in  the  upper  75  m.  We 
use  an  empirical  formula  for  ct  from  Mamayev  (1975); 


c  *  28. 152-0. 0735T-0.00469T2+(0. 802-0. 002T) (S-35)  , 


(2.1) 


Figure  2.1  Left:  Profiles  of  temperature  and 
salinity  at  49°35‘N,  145°0B’W,  0845 
GMT,  21  August  1977  (solid)  and  49° 
38'N,  145°08'W,  1320  GMT,  24  August 
1977  (dotted).  From  Davis  et  al. 
(1981)  . 

Figure  2.2  Right:  T-S  relation  from  CTD  casts 
averaged  on  19  (solid) ,  24  (dashed^ , 
and  28  (dotted)  August,  1977.  The 
6  September  curve  (not  shown)  closely 
coincides  with  that  from  19  August. 
From  Davis  et  al.  (1981). 


where  T  is  temperature  in  cC,  and  S  is  salinity  in  °/oo 
(parts  per  thousand).  We  are  interested  in  obtaining  the 
density  gradient. 


We  write 


dS  ~ 
di  = 

/dS\  dT 

\m)  Jz  '  s 

=  3  , 

(2.2) 

where  S’ 

=  32 

.  7°/oo  and 

(dS/dT)  = 

-0 . 041 ( °/oo/ 0 C )  are 

representative 

values  of 

salinity  and  salinity  gradient. 

estimated 

from 

Figure  2.2. 

The  density 

gradient  becomes 

dafc  j, 

dz 

-(0.1017+0. 

dT 

•0093T)  . 

a  z 

(2.3) 

Then,  from  the  definition  of  Vaisala  frequency 


n2  =  -  a.  |£  , 

p  dz 


we  obtain 


2  ~ 


N  = 


(9.93+0.91T) «10'4  || 


(sec 


-2, 


(2.4) 


(2.5) 


In  practice,  we  use  pairs  of  current  meters  to  compute  N^. 
With  two  current  meters  1  and  2  with  vertical  separation  Az, 
T  is  the  average  (Ti+T2)/2,  and  dT/dz  is  the  first  dif¬ 
ference  ( Tj-T2 ) /Az . 


Section  3 

TEMPORAL  VARIABILITY  OF  SHEAR  DISTRIBUTION 


In  this  section  we  examine  the  variability  of 
shear.  Much  of  the  shear  variance  is  oscillatory  at  the 
inertial  frequency.  This  is  shown  clearly  in  Figure  3.1, 
which  shows  several  examples  of  rotary  spectra  of  shear.  The 
spectral  peak  at  the  local  inertial  frequency  is  the  dominant 
spectral  feature  in  all  cases.  The  inertial  frequency  is 
f  *  2D  sin  \,  where  0  is  the  earth's  angular  rotation  rate 
and  X  is  the  latitude. 

In  Section  3.1  we  describe  the  process  of  complex 
demodulation,  and  how  it  can  be  used  to  determine  the  fre¬ 
quency  and  angle  of  propagation  of  an  internal  wave  packet. 
In  Section  3.2  we  discuss  the  bimodal  probability  density 
functions  observed  for  inertial  waves  over  short  term 
intervals.  As  we  increase  the  time  interval  of  observation, 
the  bimodal  probability  density  becomes  smeared,  and  within 
about  three  inertial  periods,  approaches  Gaussianity.  This 
result  is  shown  in  Section  3.3. 

3.1  COMPLEX  DEMODULATION 

The  process  of  complex  demodulation  of  a  shear 
time  series  can  give  us  the  temporal  variability  of  near- 
inertial  shear  amplitude  and  phase.  Let  Uz(t)  *  uz+ivz 
be  a  vector  time  series  of  shear.  We  seek  a  complex  time 
series  D(t)  such  that  D(t)exp( -ift)  most  closely  resembles 
Uz  in  a  time  interval  2T  centered  around  t,  where  T  is  the 
inertial  period  2u/f.  Perkins  (1970)  showed  that 
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torclockwi 


D(t) 


1 

2T 


(3.1  ) 


t+T 

•r  exp(  i  ft)  U  ( t )  dt  . 
t-T  '  Z 


Thus,  | D |  and  arg(D)  are  the  amplitude  and  phase  of  inertial 
shear  at  time  t.  Figure  3.2  is  an  example  of  complex  de¬ 
modulated  velocity  difference  between  29  and  32  m.  The 
dashed  curve  is  the  amplitude  and  the  solid  curve  is  the 
phase.  A  slow  decrease  in  phase  with  time  corresponds  to  a 
frequency  slightly  higher  than  inertial.  For  example,  be¬ 
tween  days  4  and  8,  the  amplitude  increases  but  slowly,  while 
the  phase  decreases  by  about  120'.  This  corresponds  to  a 
frequency  increment  iw  =  .003  cph,  giving  a  total  frequency 
f  +  iu  *  .067  cph,  about  5%  higher  than  f.  This  frequency 
increment  is  of  the  same  order  as  observed  b}  others  (e.g., 
Kundu  ,  1976 ) . 

In  order  that  internal  waves  be  able  to  propagate, 
their  frequency  to  must  be  in  the  range  f  <  w  <  N.  As  to 
decreases  toward  f,  the  group  velocity  vector  approaches  the 
horizontal  direction,  and  the  vertical  component  vanishes  at 
to  *  f.  The  angle  of  propagation  ?  of  the  group  velocity 
vector  with  respect  to  the  horizontal  is  approximately 

e  =  y/2(  w/f-1)  ,  (3.2) 

where  the  assumptions  to-f  <<  f  <<  N(z)  have  been  made  (see 
Rubenstein,  1981).  In  the  present  case  illustrated  from 
days  4  to  8  in  Figure  3.2,  the  rms  value  of  N(z)  is 
0.023  sec-1.  Therefore  the  angle  of  propagation  is 
1.6x10-3  rad,  or  about  0.1*. 
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Figure  3.2  Complex  demodulation  of  velocity  difference 

between  29  and  32  m.  Dashed  curve:  amplitude, 
solid  curve:  phase.  When  the  slope  in  phase 
•  is  positive  (negative)  the  frequency  is  lower 

(higher)  than  the  inertial  frequency.  From 
equation  (3.1)  in  the  text. 
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DISTRIBUTION  OF  INERTIAL  OSCILLATIONS 


Over  time  scales  of  order  one  inertial  period,  the 
probability  density  of  shear  is  a  bimodal  distribution,  while 
over  much  longer  time  scales,  the  probability  density  ap¬ 
proaches  that  of  a  Gaussian  distribution.  In  an  ideal  case 
with  a  constant  amplitude  inertial  oscillation  with  zero 
phase  shift,  the  shear  vector  traces  out  a  circular  pattern. 
Each  shear  component  has  a  bimodal  probability  density  func¬ 
tion.  With  the  introduction  of  a  component  of  higher  fre¬ 
quency  shear,  the  circular  pattern  becomes  smeared.  If  the 
amplitude  of  the  inertial  oscillations  varies  slowly,  then 
this  low  frequency  component  of  shear  smears  the  circular 
pattern  still  further.  If  these  "noise"  components  are  suf¬ 
ficiently  strong,  then  over  the  course  of  several  inertial 
periods,  the  probability  density  function  approaches 
Gaussianity . 

We  can  examine  this  process  by  developing  a  simple 
statistical  model  of  "noisy"  inertial  shear.  We  begin  with 
the  ideal  case,  by  stipulating  that  the  shear  vector  oscil¬ 
lates  in  a  circle,  and  each  of  the  shear  components  oscillate 
sinusoidally,  but  out  of  phase  by  ti/2.  For  example,  let 

xi  =  x^sinft  (3.3) 


be  the  uz  component  of  inertial  shear,  where  Xi  is  the 
amplitude  of  the  oscillation.  Then  it  is  easy  to  show  that 
the  probability  density  function  is 


f(x.) 


/ 


0 


xi|  >  Xi  .  (3.4) 
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Now  suppose  that  we  add  to  this  inertial  shear  an 
extra  "noise"  component  xn.  We  let  this  component  have  a 
uniform  density  function: 


l/(2Xn) 


9U„>  - 


lx  |  < 

n '  —  n 

ix  i  >  Xr,‘ 

n  i  n 


(3.5) 


We  are  interested  in  the  probability  density  h  of  the  total 
shear  x  =  x^  +  xn.  We  compute  h  by  performing  the 
following  integral; 


h( x )  =  j  f(x-xn)  g(xn)dxn> 


(3.6) 


We  combine  equations  (3.4)— (3.6)  and  after  a  bit  of  algebra 
we  obtain  a  symmetric  probability  density  function.  The 
positive  side  of  this  function  is  written  as  follows; 


1  1  /Xr,*X 

h(  x )  -.(2*  Xn)-1  sin 


-1  /Xn+X 

+  sin  ^  x' 


(3.7a) 


in  the  range  0  <.  x  £  xi  ~  Xn,  and 


*  (2nxn)-1  +  | 


(3.7b) 


in  the  range  Xi  -  xn  £  x  S.  xi  +  xn  •  Figure  3.3  shows 

h ( x )  for  Xi  =  1,  Xn  =  0.1,  0.3,  and  1.0.  For  Xn  =  0.1, 

the  probability  density  function  is  strongly  bimodal;  the  two 
cusps  are  at  x  =  +  (Xi*Xn).  As  xn  increases,  the  cusps 

approach  each  other,  and  converge  when  Xn  =  1;  that  is  to 

say,  the  distribution  becomes  unimodal. 


PROBAB  I  LI  I  Y  Di  NS  I  1 


SriEAR 


Figure  3.3  Probability  density  functions  h(x)  given  by 

equation  (3.7).  The  ratios  Xn/Xi  =  0.1,  0.3, 
and  1.0  are  proportional  to  the  ratio  of  the 
variances  of  the  noise  and  inertial  oscillation 
distributions. 
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Over  a  short  duration,  this  model  is  often  a 
reasonable  approximation  to  the  observed  probability  density. 
Figure  3.4  shows  a  trace  of  the  vector  of  velocity  difference 
between  26  and  38  m,  over  an  inertial  period  ( ~1 6  hours). 
The  trace  appears  to  be  approximately  circular,  with  seem¬ 
ingly  random  noise  superimposed.  Figure  3.5  shows  the  dis¬ 
tribution  of  the  u-component  of  this  vector  time  series.  The 
agreement  with  the  model  probability  density  function  from 
(3.7),  with  Xn  =  0.3,  appears  to  be  very  good.  We  do  not 
attempt  to  optimize  or  to  quantify  this  agreement.  We  only 
wish  to  point  out  that  the  observed  distribution  is  bimodal, 
and  can  be  understood  on  a  qualitative  level  as  the  com¬ 
bination  of  the  "oscillatory"  distribution  (3.4)  and  the 
uniform  distribution  (3.5).  In  most  of  the  cases  where  shear 
component  time  series  were  analyzed  over  an  inertial  period, 
the  distributions  were  found  to  be  bimodal.  However,  in  many 
cases  the  distributions  did  not  compare  favorably  with  (3.7), 
because  the  distributions  had  significantly  large  skewness 
levels.  In  these  cases,  non-inertial  frequency  components 
were  significant. 

3.3  LONG-TERM  GAUSSIAN  DISTRIBUTION 

Over  a  duration  of  several  inertial  periods,  low 
frequency  shear  components  accumulate  and  tend  to  smear  out 
the  bimodal  distributions.  A  duration  of  three  inertial 
periods  is  sufficient,  generally,  to  generate  nearly  Gaussian 
distributions.  We  have  analyzed  the  skewness,  kurtosis,  and 
Lilliefors  test  statistic  for  48  hr  shear  distributions  over 
Az  =  3m.  The  48  hr  interval  is  almost  exactly  three  15.7  hr 
inertial  periods.  The  shear  was  sampled  at  4  hr  intervals, 
in  order  to  achieve  statistically  independent  samples. 
Section  5.1.1  presents  the  rationale  for  this  choice.  Also, 
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Figure  3.4  Vector  trace  of  velocity  difference 

between  depths  26  and  38  m,  during  an 
inertial  period,  about  16  hours.  Over 
the  course  of  an  inertial  period,  the 
trace  is  approximately  circular,  with 
higher  frequency  fluctuations  superposed, 


3-10 


j n 


- ' - - - 1 - - - 1 - - - r - - - 1 - - -  ’ - - - r  - - 1 - - - 1 - ^ 

-i  30  -iM  -.08  —0  72  -0o6  0  00  j.jt  0?i  VO?  ‘  ^  .  ?0 

<5^/6z  (“0';  sec"') 


x  x  x  =  Theoret.cci  D .  i  tr,  Dot  r-r. 
Bars  =  Ocjer.fej 


Figure  3-5  Shear  distribution  over  depth  interval  26-38  m, 
over  an  inertial  period,  about  16  hours.  The 
theoretical  distribution  function  from  equation 
(3.7)  is  bimodal,  with  Xn  =  0.3.  The  fit  has  not 
been  optimized,  but  does  appear  to  be  reasonable. 
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we  refer  the  reader  to  Section  5.2  for  a  discussion  of  the 
Lilliefors  test  for  Gaussianity.  Skewness  is  defined  as 
j  ( x  i  -T)  3/n  ^3  /  and  kurtosis  is  defined  as  r(xi~‘x)4/ 
Na^,  where  T  =  vXj/N  and  c 2  =  v(Xi-)T)2/N  are  the 

mean  and  variance,  respectively.  The  expected  skewness  of  a 
Gaussian  distribution  (as  well  as  any  symmetric  distribution) 
is  zero.  The  expected  kurtosis  of  a  Gaussian  distribution  is 
3.  For  this  reason,  we  have  subtracted  3  from  each  kurtosis 
statistic. 


Figure  3.6  shows  an  example  of  the  analysis 
results.  This  example  has  been  selected  to  be  representative 
of  the  entire  body  of  results.  With  but  few  exceptions, 
skewness  and  kurtosis-3  oscillate  about  zero,  within  the 
limits  +1.5.  The  test  statistic  generally  falls  below  0.24, 
the  95%  confidence  (a  =  .05)  level.  The  implication  is  that 
these  short  samples  come  from  a  population  which  may  reason¬ 
ably  be  approximated  as  a  Gaussian  distribution. 
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Section  4 

STATISTICAL  MODEL  OF  MEAN-SQUARE  SHEAR 


In  this  section  we  derive  a  simple  spectral  model 
of  shear  to  predict  how  measurements  of  mean-square  shear 
vary  with  the  vertical  separation  Az  between  current  meters. 
This  analysis  is  largely  taken  from  Grabowski  (1981),  but 
the  deviation  is  slightly  different.  Then  we  compare  these 
predictions  with  measurements  in  the  seasonal  thermocline 
from  the  MILE-1  current  meters.  The  comparison  is  favorable, 
and  lends  support  to  the  composite  shear  spectrum  of  Gargett 
et  al .  (1981),  upon  which  our  model  is  based. 

4.1  PREDICTIONS  OF  MEAN-SQUARE  SHEAR 


We  begin  by  representing  the  complex  valued 
velocity  field  U(z)  =  u(z)  +  iv(z)  as  a  Fourier  integral 


-/ 


U(k)  eikz  dk  , 


(4.1) 


where  k  is  the  vertical  wavenumber.  The  one-sided  velocity 
power  spectral  density  function  (defined  over  positive  k)  is 


Pu(  k)  =  2<U( k) U  ( k) >  , 


(4.2) 


where  the  brackets  denote  an  ensemble  average. 

Suppose  we  were  to  sample  U(z)  at  discrete  inter¬ 
vals,  and  estimate  the  shear  S(z;Az)  by  performing  a  first 
differencing  operation; 


S( z ; Az )  *  U ( z  +  Az )  -  U ( z ) 
AZ 


4-1 


Taking  the  Fourier  transform  of  (4.3)  we  get 


r 


*- 


» 


» 


t 


♦ 


t 


S(k;Az)  =  -2iU(k)  exp(-ikAz/2)  sin(kAz/2)  .  (4.4) 

A  z 

Therefore,  the  power  spectral  density  of  shear  Ps(k;Az) 
over  a  finite  interval  Az  is 

Pg(k;Az)  =  4Pu(k)  ( sin  ( 


kAz/2 )\ 2  .  (4.5) 

Az  / 


In  the  limit  as  Az  goes  to  zero, 


P  (k)  =  lim  P  (k;A  z)  *  k‘p  (k) 
s  Az-0  S  u 


(4.6) 


In  Figure  4.1  we  illustrate  a  shear  spectrum  model 
Ps(k),  after  the  composite  spectrum  by  Gargett  et  al.  The 
form  of  this  model  is 


PsOO 


(NX 

|N2/k 


k  <  ko 
kQ  <  k  <  kb. 


(4.7) 


where  N  is  the  local  Vaisala  frequency,  kQ  is  taken  as 
0.1  cpm,  and  the  buoyancy  wavenumber,  k^,  is  taken  as 
1  cpm.  we  do  not  include  contributions  from  the  high 
wavenumber  turbulence  dissipation  range  in  this  model.  We 
substitute  (4.6)  and  (4.7)  into  (4.5)  and  integrate  over 
wavenumber  to  obtain  the  expected  value  of  mean  square 
shear; 
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k. 
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Vertical  Wavenumber  k  (cpm) 


4.1  Idealized  shear  spectrum  after  Gargett 
et  al.  (1981).  The  rolloff  wavenumber 
is  kQ  =  .1  cpm,  and  the  buoyancy  wave- 
number  is  kh  =  1  cpm. 


<  S( z ; Az)  >  =  N  /k  J  sin(kAz/2) 

o  „  — r  .  _  To - 


k  Az/2 


+  n2  /  k_i n 

ko  L 


sin(kAz/2)  dk. 

■  Tun — 


Writing  y  =  kAz/2,  y0  =  k0Az/2,  and  y^  =  kbAz/2 

(4.8)  becomes 


<|s(z; Az) J2>  = 


Us  /°  si°2y  dy  +  n2  /b  sin2y  dy 

°  o  y2  yo  y3 


Grabowski  (1981)  evaluated  the  first  integral  1^  in  (4.9); 


I,  -  COS  2yo  "  1  +  Si ( 2y_ ) 

2jT 


(4. 1C 


and  the  second  integral  I2  is 


■*[ 


y~2  (cos2yb  -  1)  -  2y~2  sin2yt 


(cos2y0  -  1)  +  ay;1  sin2yJJ 


where 


-Ci(2y0)  +  Ci(2yb)  , 


(4.11 


*  —1 

Si(x)  »  /  t  sint  dt 
0 


(4.12) 


Cl  ( x )  =  -J 
x 


cost  dt 


are  the  sine  and  cosine  integrals.  Model  predictions  of 
normalized  mean-square  velocity  differences 


<  S(zjAz)  2>  Az2k 


o  =  Az(Ix  +  kQAzI2/2) 


(4.13) 


are  tabulated  in  Table  4.1.  From  Table  4.1  it  is  evident 
that  mean  square  velocity  differences  become  asymptotically 
proportional  to  Az  for  large  Az.  in  the  next  section,  we 
present  a  graph  of  the  relationship  (4.13),  and  compare  with 
the  values  of  mean-square  velocity  differences  from  the  MILE 
current  data. 

Table  4.1 

Mean-Square  Velocity  Differences 


Az  (m) 


<  S(  z;  Az)  2>  Az2  k0/2N 2 


2.15 

4.84 

8.17 

11.39 
14.32 
17.24 

20.40 
23.70 


26.89 


W  «af>r' 
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STATISTICAL  MODEL-DATA  COMPARISONS 


In  th  s  section  we  compare  mean-square  shear  with 
our  model  prediction.  We  begin  by  showing  that  S2,  computed 
from  current  meters  over  each  of  four  different  vertical 
separations  Az,  is  linearly  proportional  to  N2.  The  overbars 
denote  temporal  averages  over  a  19-day  period.  Using  the 
method  of  least-squares,  we  compute  estimates  of  the 
normalized  mean-square  velocity  difference  kQ( S2/2N 2) Az2 , 
as  a  function  of  Az.  Finally  we  compare  these  results  with 
the  statistical  model  prediction  (4.13). 

Figures  4. 2-4. 5  show  observed  values  of  S2  as 
functions  of  N2,  for  Az  =  3,  6,  12,  and  21  m.  Mean-square 
total  shear  S2  is  denoted  by  triangles,  and  is  plotted 
against  the  left-hand  axes.  The  individual  components  of 
mean-square  shear,  ( Au/  Az ) 2  and  (  Av/ Az ) 2 ,  are  plotted  on 
twice  the  scale  as  the  total  mean-square  shear.  The  dashed 
lines  show  the  origin-forced  least-squares  best  fits;  S2  * 
bN2.  In  each  figure,  the  correlation  coefficient  is  better 
than  0.99. 


Table  4.2  lists  the  coefficients  b  of  the  least- 

squares  regressions  for  each  value  of  Az.  The  value  of  b 

decreases  with  Az ,  and  is  approximately  proportional  to 

Az“l .  That  is  to  say,  mean-square  shear  becomes 

proportional  to  A z_i,  and  mean-square  velocity 
differences  S2az2  become  proportional  to  Az,  for  large 
Az . 
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Figure  4 


MEAN  S**2  VS.  N**2.  DZ  =  3  M 


MEAN  N*»2  (E— a  SEC»»— 2) 


TRIANGLES:  TOTAL  s2 

_ 

CROSSES:  2x  (DU/DZ)  ,  (DV/DZ) 

DASHED  LINE:  S2  =  0.6IN2 


.2  F  vs.  N^  from  eight  pairs  of  current  meters 
in  the  depth  range  26-50  m,  with  -z  -  3  m. 
Individual  shear  components  (crosses)  are  plotted 
on  twice  the  scale  of  the  total  S^  (triangles). 
Dashed  line  is  best  fit  through  the  triangles. 
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Figure 
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p  0  2.0  4.0  6  0  8  0  10  0 

MEAN  N«»2  (E-a  SEC»»— 2; 


TRIANGLES:  TOTAL  S2 

CROSSES:  2x  (DU/DZ)  ,  (DV/DZ)2 

DASHED  LINE:  S2  =  0.46N2 


4.3  S2  vs.  N2  from  seven  pairs  of  current  meters 
in  the  depth  range  26-50  m,  with  Lz  -  6  m. 
Individual  shear  components  ( crosses )  are  plotted 
on  twice  the  scale  of  the  total  S*  (triangles). 
Dashed  line  is  best  fit  through  the  triangles. 
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MEANi  N»*2  SEC*»-2) 


2 

TRIANGLES:  TOTAL  S 

— - j  -  j 

CROSSES:  2x  (DU/DZ)  .  (DV/DZ) 

DASHED  LINE:  S2  =  0.29N2 


Figure  4.4  S^  vs.  N^  from  five  pairs  of  current  meters 
in  the  depth  range  26-50  m,  with  L z  =  12  m. 
Individual  shear  components  (crosses)  are  plotted 
on  twice  the  scale  of  the  total  S^  (triangles). 
Dashed  line  is  best  fit  through  the  triangles. 
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Figure  4.5  S^  vs.  N^  from  three  pairs  of  current  meters 
in  the  depth  range  26-70  m,  with  Lz  =  21  nt. 
Individual  shear  components  (crosses)  are  plotted 
on  twice  the  scale  of  the  total  S^  (triangles). 
Dashed  line  is  best  fit  through  the  triangles. 


Table  4.2 

Least-Squares  Regression  Coefficients  =  bN^ 


Az  (m) 

b  =  S2/N2 

o  T  T  2 

b  Az^  =  (S  /Nz)  Az 

k  — 2 
*OS^  AZ 

2 

2N 

3 

0.61 

5.49 

1.72 

6 

0.46 

16.56 

5.20 

12 

0.29 

41.76 

13.12 

21 

0.16 

70.56 

22.17 

This  result  is  analogous  to  that  predicted  by  the 
statistical  model,  eq.  (4.13).  Figure  4.6  shows  the 
statistical  model  of  mean-square  velocity  differences  (solid 
curve),  the  model  times  the  factors  2±W2  (dashed),  and 
observed  values  (triangles).  The  error  bars  indicate  +2 
standard  deviations,  and  represent  the  uncertainties  in  the 
coefficients  b.  The  model  prediction  is  greater  than  the 
observed  values.  Both  the  model  and  the  observed  values  of 
mean-square  velocity  difference  increase  approximately 
linearly  with  Az,  with  slopes  1.5  and  1.16,  respectively. 

Several  possible  explanations  may  be  advanced  for 
this  discrepancy.  The  model  may  be  overestimating  the  low 
wavenumber  (k  £  ko)  spectral  level,  N2/k0.  As  the  energy 
level  of  internal  waves  is  universal  within  a  factor  of  two 
(Nunk,  1982),  and  much  of  the  shear  variance  is  due  to  low 
frequency  internal  waves,  it  is  reasonable  to  expect  shear 
spectral  level  to  vary  by  a  similar  factor.  The  dashed 
curves  in  Fig.  4.6  show  the  statistical  model  predictions 
with  the  shear  level  adjusted  by  the  factors  2±1/2. 
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A  =  OBSERVED 

-  =  MODEL 

-  =  (2±1/2)  MODEL 

♦  '  ' 

Figure  4.6  Normalized  mean-square  velocity  difference,  as  a 
function  of  vertical  separation.  The  obser¬ 
vations  are  based  on  MILE  current  meter  pairs  in 
the  seasonal  thermocline,  in  the  depth  range  26- 
9  70  m.  The  dashed  curves  show  the  model  adjusted 

by  the  factors  2*1/2. 
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The  assumed  shape  of  the  shear  spectrum  may  not  be 
exact.  Gargett  et  al .  (1981)  estimate  the  spectral  roll-off 
wavenumber  kD  to  be  0.1  cpm,  and  claim  this  number  to  be 
correct  within  a  factor  of  two.  However,  our  statistical 
model  is  not  sensitive  to  changes  in  kQ,  and  varying  kQ 
by  a  factor  of  two  will  not  produce  better  agreement  with  the 
observations.  Gargett  et  al .  mention  that  the  low  wavenumber 
portion  of  the  shear  spectrum  might  not  be  truly  flat,  and 
that  Sanford  suggests  a  change  in  spectral  slope  from  k®  to 
k“*5  around  k=.02-.03  cpm.  If  this  spectral  shape  were 
correct,  then  the  low-wavenumber  portion  of  the  spectrum 
would  gain  in  relative  importance.  As  a  result,  the  model's 
prediction  of  mean-square  shear  would  increase  prefer¬ 
entially  for  large  values  of  A  z. 

Another  possibility  is  that  N2  is  not  computed 
correctly.  We  assumed  in  Section  2.2  that  fluctuations  in 
salinity  gradient  are  nearly  proportional  to  fluctuations  in 
temperature  gradient,  and  that  the  salinity  gradient's 
contribution  to  the  total  density  gradient  is  small.  If  this 
assumption  is  not  valid,  then  our  estimate  for  may  be 

faulty.  For  example,  if  the  fluctuations  in  salinity 

gradient  were  less  extreme  than  those  in  temperature 
gradient,  then  we  have  overestimated  N2. 

Yet  another  possible  explanation  for  the  dis¬ 
crepancy  lies  in  the  fact  that  N2  is  a  strongly  varying 

function  of  depth  z.  When  Az  is  as  large  as  21  m,  the 

local  values  of  N2  at  depths  z  and  z+Az  may  differ 

significantly  from  the  average  value  over  this  depth 
range.  This  may  cause  a  problem  since  S2  is  normalized  by 
the  average  value  of  N2. 
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Our  statistical  model  provides  a  framework  for 
obtaining  information  about  the  vertical  length  scales  of 
shear  from  current  meter  time  series.  There  is  good  qualita¬ 
tive  agreement  between  values  of  shear  computed  from  current 
meter  time  series  and  those  computed  from  the  statistical 
model.  The  model  is  derived  from  a  vertical  wavenumber 
spectrum  which,  in  turn,  is  based  on  vertical  profiles. 
This  is  important  because  it  is  difficult  to  compute  reliable 
vertical  wavenumber  shear  spectra  from  vertical  profiles  in 
the  statistically  nonstationary  regime  of  the  seasonal 
thermocline . 


Section  5 

PROBABILITY  DISTRIBUTION  OF  SHEAR 
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In  this  section  we  show  that  shear  distribution, 
over  a  time  interval  long  compared  to  an  inertial  period, 
may  be  reasonably  approximated  by  a  Gaussian  function.  in 
Section  5.1  we  show  that  shear  time  series  should  be  sampled 
at  4  hr  intervals  (1/4  inertial  period)  in  order  to  obtain 
statistically  independent  (uncorrelated)  samples.  We  also 
show,  using  the  Smirnov  test,  that  the  samples  of  Au/Az  and 
Av/Az  shear  components  come  from  parent  populations  whose 
distributions  are  identical. 

in  Section  5.2  we  use  the  Lilliefors  test  to  show 
that  Au/Az  and  Av/Az  have  distributions  that  are  approxi¬ 
mately  Gaussian.  Then  in  Section  5.3  we  apply  the  Kolmogorov 
test.  We  show  that  the  Au/Az  and  Av/Az  distributions  agree 
well  with  the  particular  Gaussian  distributions  whose 
variance  parameters  are  based  on  the  statistical  model  for 
S^(N^", Az)  from  Section  4. 

5.1  PRELIMINARY  STATISTICAL  TESTS 

5.1.1  Autocorrelation  and  Run  Test  for  Randomness 


In  order  to  perform  quantitative  statistical  tests 
on  the  shear  probability  distribution,  we  must  first  realize 
and  account  for  the  fact  that  shear  measurements  are  not 
independent,  but  are  autocorrelated .  We  use  two  methods  to 
determine  autocorrelation  lag  interval.  The  first  method 
involves  estimating  this  interval  directly  from  shear  auto¬ 
correlation  functions,  and  the  second  involves  the  use  of  the 
run  test. 


5-1 


In  general,  shear  in  the  shallowest  region  of  the 
seasonal  thermocline  rapidly  becomes  uncorrelated,  while 
shear  in  the  slightly  deeper  regions  remains  correlated  for  a 
length  of  time  of  order  one  inertial  period.  Figures  5.1  and 
5.2  are  two  examples  of  autocorrelation  functions  of  shear 
between  current  meter  pairs  at  depths  29-32  m  and  41-44  m, 
respectively.  In  Figure  5.1,  the  shear  quickly  becomes 
uncorrelated.  After  a  lag  of  1  hour,  the  shear  reaches  half¬ 
correlation.  At  larger  lags,  the  function  continues  to 
oscillate  at  the  inertial  frequency.  The  amplitude  of  these 
oscillations  is  about  +  0.2.  Figure  5.2  shows  a  somewhat 
different  behavior.  The  inertial  frequency  oscillations 
decay  more  slowly,  and  reach  half-correlation  after  a  lag  of 
about  20  hours.  These  examples  have  been  chosen  to  be 
representative  of  the  entire  body  of  results. 

In  all  of  the  autocorrelation  results,  the  shear 
autocorrelation  function  first  passed  through  zero  at  a  lag 
of  about  4  hours.  This  result  is  obvious,*  since  inertial 
oscillations  are  an  important  component  of  the  shear,  the 
signal  reaches  zero  correlation  after  1/4  of  an  inertial 
period  (15.66  hours),  or  about  4  hours. 

A  second  test  was  also  used  to  ascertain  the  degree 
of  randomness;  the  run  test.  We  individually  apply  this  test 
to  each  of  the  horizontal  components  of  vertical  shear.  We 
begin  this  test  by  subtracting  the  mean  values  from  both 
shear  components.  A  run  is  defined  as  a  sequence  of  either 
positive  or  negative-valued  shear  observations  that  is  fol¬ 
lowed  by  observations  of  the  opposite  sign  (or  no  observation 
at  all) .  The  total  number  of  runs  in  a  time  series  gives  us 
an  indication  of  its  randomness.  We  searched  for  the  minimum 
subsampling  interval  for  which  the  run  test  hypothesis  could 
be  accepted,  consistently ,  at  the  a  =  0.20  significance 
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Figure  5.1  Autocorrelation  of  shear  over  depth  interval 

29-32  in.  The  solid  cur'.'e  is  the  autocorrelation 
of  du / Lz  ,  and  the  dashed  curve  is  for  Lv/Lz. 
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Figure  5.2  Autocorrelation  of  shear  over  depth  interval 

41-44  m.  The  solid  curve  is  the  autocorrelation 
of  Au/Az,  and  the  dashed  curve  is  for  Lv/lz. 


level.  At  this  level,  sample  populations  of  a  truly  random 
parent  population  should  fail  the  run  test  less  than  20%  of 
the  time.  (A  detailed  description  of  the  run  test  may  be 
found  in  Conover,  1971.)  We  find  that  the  minimum  sub¬ 
sampling  interval  is  four  hours. 
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This  means  that  if  we  subsample  at  the  rate  of  once 
per  four  hours,  adjacent  pairs  of  observations  are  nearly 
equally  likely  to  be  of  the  opposite  and  of  the  same  signs. 
This  result  is  in  agreement  with  the  autocorrelation  func¬ 
tion,  which  first  passes  through  zero  at  4  hours.  We  choose 
to  use  this  definition  of  randomness,  or  "statistical  inde¬ 
pendence."  In  all  of  the  statistical  tests  in  the  remainder 
of  Section  5,  we  subsample  at  4  hour  intervals.  If  the 
inertial  amplitude  remains  correlated  over  a  significantly 
longer  interval,  then  all  of  the  samples  will  not  be  statis¬ 
tically  independent,  and  the  number  of  independent  samples 
will  be  overestimated.  As  a  result,  the  tendency  will  be  for 
a  true  null  hypothesis  to  be  rejected. 

5.1.2  Smirnov  Test  for  Identical  Distributions 

Before  proceeding  with  the  tests  for  Gaussianity, 
we  first  determine  whether  or  not  the  distribution  functions 
associated  with  the  two  populations  of  Au/Az  and  Av/Az 
observations  are  identical.  For  this  purpose  we  use  the  two- 
sided  Smirnov  test  described  by  Conover  (1971).  The  test 
statistic  T  is  defined  to  be 

T  =  s  ijp  j  S 1  ( x )  -  S2(x)  |  ,  (5.1) 

where  Si(x)  and  S2(x)  are  the  observed  distribution 
functions  of  Au/Az  and  Av/Az,  respectively.  The  symbol  sup 


» 
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denotes  supremum,  or  maximum  value  over  the  range  of  x.  In 
other  words,  the  test  statistic  is  the  greatest  vertical 
distance  between  the  distribution  functions  Si  ,x)  and 
S2(x)  . 

We  analyzed  all  of  the  shear  distributions  in  the 
upper  thermocline,  Au/Az  and  Av/Az,  for  Az  =  3,  6,  12,  and 
21  m.  We  first  subsampled  at  4  hour  intervals,  and  then  sub¬ 
tracted  the  means  computed  over  the  entire  19  day  duration  of 
the  time  series.  The  two  resulting  time  series,  each  approx¬ 
imately  of  length  19x24/4  =  114  points,  were  sorted  and  the 
distribution  functions  computed.  The  Smirnov  test  is  appli¬ 
cable  even  if  the  two  sample  populations  are  not  equal  in 
size  (but  in  our  case  the  populations  were  equal  in  size.'. 

In  all  of  the  cases  we  studied,  the  null  hypoth- 
esis--that  the  two  distributions  were  identical--was  accepted 
at  the  a  =  0.20  level  of  significance.  This  result  is  not 
very  surprising  in  light  of  Figures  4. 2-4. 5.  In  this  set  of 
figures,  the  individual  components  of  mean-square  shear, 
i.e.,  the  variances  of  Au/Az  and  Av/Az,  are  nearly  equal. 
Therefore  we  do  not  expect  significant  differences  between 
the  statistical  tests  for  au/az  and  Av/Az  in  Sections  5.2- 
5.4. 


5.2  LILLIEFORS  TEST  FOR  GAUSSIANITY 

In  this  section  we  apply  the  Lilliefors  test  for 
Gaussianity  (normality)  to  the  horizontal  components  of 
vertical  shear,  Au/Az  and  av/az.  The  null  hypothesis  is  that 
the  random  sample  has  a  Gaussian  distribution,  with  un¬ 
specified  mean  and  variance.  In  order  to  assure  randomness, 
we  subsample  the  shear  components  at  4  hour  intervals,  as 
discussed  in  Section  5.1.1. 
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Let  Xi,  i  =  1,2, ...n  be  the  measured  samples  of 
shear,  Au./  Az  or  Av/Az.  We  normalize  the  sample  values  by 
subtracting  out  the  mean  "x  and  dividing  by  the  standard 
deviation  c; 


*i  -  ILL! 

c 


(5.2) 


IE 


n  1/2 

J_  V*  (x.  -  x)z 

n-1  ZmJ 


The  new  samples  z^ ,  i=l,2,...,n  have  zero  mean  and  unity 

standard  deviation.  The  "standard"  C?=0,  oz=l)  Gaussian 
probability  density  function  is 


f(z)  =  (  2n)'1/2  exp(  -z2/2  )  , 


(5.5) 


and  the  associated  "standard"  distribution  function  is 


*  1 

F(z)  =  /  f(z)  dz  =  ■*■  +  erf(z)  , 


where  erf  is  the  error  function. 


Figure  5.3  is  a  representative  histogram,  showing 
the  observed  probability  density  of  Au/Az  computed  from  the 
current  meter  pair  at  depths  41-44  m.  The  abscissa  has  been 
normalized  as  in  equations  (5.2)-(5.4).  The  observed 
probability  density  is  compared  with  the  Gaussian  function 
(5.5),  denoted  by  x's.  The  agreement  is  appealing  to  the 
eye,  but  we  need  an  objective  test  of  the  goodness  of  fit. 
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•  Figure  5.3  Frequency  distribution  histogram  of  shear 
component  du/dz  over  41-44  m.  The  19-aay 
record  has  been  subsampled  at  4  hour 
intervals,  and  normalized  by  the  standard 
deviation.  The  crosses  denote  a  normalized 
Gaussian  function,  from  equation  (5.5). 
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Some  authors  use  the  chi-square  test  for  goodness  of  fit,  but 
the  results  depend  somewhat  or.  the  number  of  bins  chosen  for 
the  test.  Instead,  we  prefer  to  use  a  test  of  the 

"Kolmogorov-Smirnov  type",  that  .s,  one  based  on  the  dis¬ 

tribution  function.  The  Lilliefors  test  for  Gaussianity  is 
of  this  type. 

Figure  5.4  is  the  observed  probability  distribution 
of  Au/Az,  corresponding  to  the  same  probability  density- 
histogram  shown  in  Figure  5.3.  The  dashed  curve  is  the 

standard  function  distribution  (5.6).  The  Lilliefors  test 
statistic  T  is  computed  from  the  standard  distribution 

function  F(z)  in  (5.6),  and  from  the  observed  distribution 
S(z)  of  the  normalized  shear  z; 

T  =  supj  F( z )  -  S(z)  |  .  (5.7) 


Thus,  T  is  the  maximum  vertical  distance  between  the  standard 
distribution  function  (dashed  curve)  and  the  observed  dis¬ 
tribution  (solid  curve)  in  Figure  5.4.  We  can  estimate  the 
test  statistic  T  graphically.  We  find  that  for  this  case,  T 
=  0.049,  which  occurs  at  z  =  1.0.  This  distance  falls  well 
within  the  80%  confidence  interval  0.068.  Therefore  the  null 
hypothesis  (that  the  distribution  is  Gaussian)  is  accepted  at 
the  a  =  0.20  level  of  significance.  Conover  (1971)  says  that 
acceptance  of  the  null  hypothesis  does  not  necessarily  mean 
that  the  parent  population  is  Gaussian.  But  it  does  mean 
that  the  Gaussian  distribution  is  not  an  unreasonable  approx¬ 
imation  to  the  true  unknown  distribution. 

Since  we  need  to  analyze  a  large  number  of  cases, 
in  practice  we  do  not  evaluate  T  graphically,  but  we  program 
a  computer  code  to  calculate  it.  Table  5.1  lists  the 
critical  levels  of  significance  a  for  all  of  the  seasonal 
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Figure  5.4  Cumulative  distribution  of  shear  component  du/dz 
over  41-44  m.  The  19-day  record  has  been  sub¬ 
sampled  at  4  hour  intervals,  and  normalized  by 
the  standard  deviation.  The  dashed  line  is  the 
error  function  cumulative  distribution,  from 
equation  (5.6).  The  Lilliefors  test  statistic, 
defined  as  the  maximum  vertical  distance  between 
the  two  distribution  curves,  is  0.047.  This 
value  falls  well  within  the  801  confidence 
interval . 
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thermocline  shear  distributions  in  the  depth  range  26-92  m. 
The  critical  level  £  is  defined  (Conover,  1971)  as  the 
smallest  significance  level  at  which  the  null  hypothesis 
would  be  rejected  for  the  given  observation.  The  results 
show  that  the  Gaussian  distribution  is  a  reasonable 
approximation  to  the  true  distribution.  About  85%  of  the  48 
cases  were  accepted  by  the  Lilliefors  test  at  the  0.20  level 
of  significance. 


Table  5.1 

Lilliefors  Gaussianity  Test  Results 

Distribution  of  critical  levels  £  for  Au/ Az  and  Av/az  for 
3,  6,  12,  and  21  m  separations. 


AZ  = 

3  m 

6  m 

12  m 

21  m 

TOTAL 

&  < 

0.01 

2 

0 

0 

0 

2 

0.05 

>  a  > 

0.01 

0 

2 

1 

0 

3 

0.10 

>  s  > 

0.05 

0 

1 

0 

0 

1 

0.20 

>  a  > 

0.10 

0 

0 

0 

1 

1 

a  > 

0.20 

14 

11 

9 

7 

41 

TOTAL 

16 

14 

10 

8 

48 
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STATISTICAL  MODEL  APPLICATION:  THE  KOLMOGOROV  TEST 


In  Section  5.2  we  applied  the  Lilliefors  test  to 
our  shear  data,  and  showed  that  Au/Az  and  Av/Az  are  well 
described  by  Gaussian  distributions.  The  mean  and  variance 
were  not  previously  specified,  but  instead  were  estimated 
from  the  data. 
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In  this  section  we  apply  the  statistical  model 
developed  in  Section  4  to  estimate  the  variance  of  iu/iz  and 
Av/iz.  The  variance  is  estimated  from  the  least-squares  fit 
in  Figure  4.6,  according  to  the  formula 


.  Is2 

-  2b 


,2  1.851Z  -  2.8 


(5.8) 


As  we  are  only  interested  in  near-inertial  frequency 
fluctuations,  and  not  in  long-term  trends,  we  have  subtracted 
the  mean  from  each  shear  record. 

We  applied  the  Kolmogorov  goodness-o f-f it  test  to 
each  of  the  shear  records  in  the  seasonal  thermocline  in  the 
depth  range  26-70  m.  The  shear  records  over  70-92  m  were  not 
included  in  this  test,  because  the  salinity  gradient's 
contribution  to  N2  in  this  interval  is  of  comparable 
importance  to  that  of  the  temperature  gradient.  The  test 
statistic  is  computed  in  the  same  way  as  for  the  Lilliefors 
test,  eq.  (5.7).  The  difference  is  that  the  variance  is 
estimated  by  (5.8),  and  the  critical  level  values  are 
evaluated  from  a  different  set  of  tables. 

Table  5.2  lists  the  distribution  of  critical  levels 
of  significance  a  for  all  the  seasonal  thermocline  shear 
records.  In  43  out  of  46,  or  93%  of  all  cases,  the  null 
hypothesis  that  the  shear  distribution  is  Gaussian  with  the 
variance  given  by  (5.8),  is  accepted  at  the  0.20  level  of 
significance . 
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Table  5.2 

Kolmogorov  Goodness-of-Fit  Test  Results 

Distribution  of  critical  levels  a  for  au/ Az  and  Av/Az  for 
3,  6,  12,  and  21  m  separations. 


Az  =  3  m  6m  12  m  21  m  TOTAL 


- ■  - 

a  <  0.01 

0 

0 

0 

0 

0 

0.05 

>  1  >  0.01 

0 

0 

0 

0 

0 

0.10 

>  a  >  0.05 

1 

1 

0 

0 

2 

0.20 

>  a  >  0.10 

1 

0 

0 

0 

1 

z  >  0.20 

14 

13 

10 

6 

43 

TOTAL 

16 

14 

10 

6 

46 

INSTANTANEOUS  PREDICTIONS  OF  SHEAR 


In  this  Section  we  address  the  following  question: 

If  we  have  a  15-minute-average  observation 
of  local  density  stratification,  with  what 
degree  of  certitude  can  we  predict  the 
value  of  shear? 

We  have  already  shown,  in  Section  5,  that  S2  and  N2  are 
well  correlated,  statistically,  over  long  time  scales.  Over 
a  long  period  however,  N2  can  change  significantly.  There¬ 
fore  we  ask  if  it  is  possible  to  predict  S2  (over  some 
vertical  separation  Az)  from  a  short-period  average  N2 , 
say,  as  short  as  15  minutes. 

Figure  6.1  shows  a  time  series  of  15-minute  samples 
of  N2  and  s2,  for  three  days,  measured  between  depths  29 
and  32  m.  We  see  that  N2  and  S2  are  often  well  corre¬ 
lated,  even  over  short  time  scales. 

This  high  degree  of  correlation  can  be  understood, 
at  least  in  part,  in  terms  of  episodes  of  shear  instability. 
From  measurements  of  shear  and  stratification  over  a  7  m 
vertical  separation  using  the  microscale  sensing  array, 
Eriksen  (1978)  argues  that  the  flow  is  dynamically  limited  by 
shear  stability.  He  shofcs  that  the  distribution  of 
Richardson  number  Ri  =  N2/S2  falls  off  sharply  below 

Ri  =  0.25.  This  cutoff  is  evocative  of  theoretical  predic¬ 
tions  of  shear  instability.  A  likely  candidate  mechanism  is 
the  Kelvin-Helmholtz  instability,  first  observed  by  woods 
(1968)  in  the  seasonal  thermocline.  The  theoretical 


pp  of  l'i-minuto  avoraqotl 
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stability  analysis  by  Miles  and  Howard  (1964)  shows  that  the 
most  unstable  disturbance  has  a  wavelength 


X  =  7.5  h,  '6.1! 

where  h  is  the  thickness  of  a  shear  layer.  Woods  (1968) 
observed  a  sheet  and  layer  structure  in  the  temperature  and 
velocity  profiles.  He  found  that  in  the  vast  majority  of 
cases,  the  layer  thickness  h  fell  in  the  range  8<h<15  cm,  and 
that  in  all  measured  cases,  the  relationship  between  X  and  h 
fell  within  20%  of  equation  (6.1).  Thus  it  seems  plausible 
that  shear  instability  may  be  occurring  over  a  small  vertical 
scale. 


Other  mixing  mechanisms,  such  as  the  convective 
instability  proposed  by  Orlanski  and  Bryan  (1969),  may  also 
be  contributing  to  the  correlation  between  N2  and  S2 . 

This  remains  a  speculative  area;  more  observational  and  theo¬ 
retical  work  is  needed. 

For  now  though,  we  seek  only  to  quantify  the  cor¬ 
relation  between  N2  and  S2  over  a  15  minute  time  scale. 

We  quantified  the  correlation  in  two  ways.  First,  we  com¬ 
puted  the  linear  correlation  coefficient  between  N2  and 

S2  for  all  pairs  of  current  meters  in  the  seasonal  thermo- 
cline,  over  the  depth  range  between  26  and  70  m,  using 

separations  of  3,  6,  12  and  21  m.  Table  6.1  shows  the  upper 

and  lower  limits  of  the  correlation  coefficient  for  each 
vertical  separation.  For  the  most  part,  the  correlation 
coefficients  fall  in  the  range  from  0.30  to  0.54.  values 

greater  than  0.20  are  correlated  at  the  5%  level  of  sig¬ 

nificance.  Therefore  we  can  say  that  for  most  of  the  data 
records,  the  N2-S2  correlation  is  statistically  sig¬ 
nificant. 
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Table  6.1 


LOWER  AND  UPPER  LIMITS  OF  LINEAR  CORRELATION 
COEFFICIENTS  BETWEEN  15-MINUTE  AVERAGED  OBSERVATIONS 

OF  N2  AND  S2 


Correlation 


Az 

Lower  Limit 

Upper  Lim 

3  m 

.34 

.54 

6  m 

.40 

.53 

12  m 

.30 

.52 

21  m 

.15 

.46 

The  second  way  in  which  we  quantified  the  cor¬ 
relation  was  by  means  of  the  distribution  of  S2/N2  mea¬ 
surements.  Figure  6.2  shows  a  histogram  of  samples  of 
S2/N2  measurements  over  the  depth  interval  41-44  m,  and 
Figure  6.3  shows  the  corresponding  cumulative  distribution 
function.  From  Table  4.2,  the  expected  value  of  S^/N2 
over  a  3  m  interval  is  0.61.  We  can  use  Figure  6.3  to  show 
that  an  estimate  of  S2  =  0.61N2  is  accurate  within  a 
factor  of  two,  47%  of  the  time,  and  within  a  factor  of  four, 
74%  of  the  time. 

Table  6.2  summarizes  the  results  for  all  pairs  of 
current  meters  in  the  seasonal  thermocline,  between  depths  26 
and  70  m,  separated  by  3,  6,  12  and  21  m.  As  Az  increases, 
the  accuracy  with  which  we  can  predict  S2  also  increases. 
The  reason  for  this  trend  is  that  the  temporal  fluctuations 
of  N2  and  S2  tend  toward  smaller  values  relative  to 
N2  and  si,  as  Az  increases. 
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Table  6.2 

AVERAGE  PERCENTAGE  OF  15-MINUTE  N2  OBSERVATIONS 
FROM  WHICH  WE  CAN  PREDICT  S2  WITHIN  A  FACTOR  OF  2  AND  4 


Accuracy  Factor 


Az 

2 

4 

3  m 

45% 

76% 

6  m 

51% 

79% 

12  m 

54% 

79% 

21  m 

60% 

85% 
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Section  7 

SUMMARY  AND  RECOMMENDATIONS 


Our  most  important  result  is  the  validation  of  a 
statistical  model  of  mean-square  shear  first  proposed 
by  Grabowski  (1981).  The  model  is  derived  from  an  idealized 
version  of  the  vertical  wavenumber  shear  spectrum  of  Gargett 
et  al.  (1981).  This  composite  spectrum  was  assimilated  from 
vertical  profiles  of  shear  in  the  main  thermocline.  The 
model  predicts  that  S2  is  proportional  to  N^/az,  where  A z 
is  the  vertical  separation  over  which  shear  is  measured.  We 
compare  the  model's  predictions  for  mean-square  velocity 
differences  with  differences  calculated  from  various  pairs  of 
MILE  current  meters,  and  find  qualitative  agreement. 

In  general,  S2  is  proportional  to  N2 ,  but  in 
the  seasonal  thermocline,  both  quantities  vary  rapidly  with 
depth.  Therefore,  vertical  shear  profilers  cannot  measure 
statistically  reliable  spectra  in  the  seasonal  thermocline. 
However,  our  statistical  model  allows  us  to  relate  shear 
spectra  obtained  from  deep  vertical  profiles  to  current 
meter-derived  S^,  as  a  function  of  and  Az.  Thus  we 
can  extend  the  range  of  validity  of  the  Gargett  et  al.  (1981) 
composite  spectrum  into  the  seasonal  thermocline. 

Inertial  oscillations  tend  to  dominate  the  higher- 
frequency  fluctuations  of  the  shear  vector.  Therefore,  over 
the  course  of  an  inertial  period,  the  probability  dis¬ 
tributions  of  Au/Az  and  Av/Az  frequently  are  bimodal.  Over 
the  course  of  several  inertial  oscillations,  slow  variations 
in  the  inertial  amplitude  tend  to  smear  out  the  probability 
distributions,  and  they  approach  Gaussianity. 


We  have  shown  that  mean  square  shear  is  nearly 
isotropic.  Figures  4. 2-4. 6  show  that  (  *  u/  Lz )  2  *{ "'xr.7  Lz )  2  ® 
S2/2 ,  where  the  overbars  denote  19  day  averages.  We  used 
our  statistical  model  for  "s2  to  parameterize  the  Gaussian 
distributions  of  the  individual  shear  components.  These 
parameterized  distributions  agree  very  well  with  the  observed 
d istributions. 

Shear  and  stratification  are  significantly 
correlated  over  short  time  scales,  as  well.  Fifteen-minute 
averaged  samples  of  N2  and  S2  measured  over  iz=3  m  have 
correlation  coefficients  in  the  range  0.34-0.54.  We  showed 
that  knowing  N2  over  a  3  m  vertical  interval  allows  us  to 
estimate  S2  within  a  factor  of  two,  46%  of  the  time,  and 
within  a  factor  of  four,  75%  of  the  time.  Similar  results 
were  obtained  over  wider  vertical  intervals. 

It  is  not  clear  from  the  MILE  current  meter  records 
whether  shear  variability  is  primarily  due  to 

•  temporal  variability, 

•  vertical  variability  due  to  vertical  dis¬ 

placements  of  the  current  meters,  or 

/ 

V 

\ 

•  vertical  variability  due  to  vertical  dis¬ 

placements  of  the  shear  field  by  internal  waves 
and  Ekman  pumping, 

• 

or  some  combination  of  these  factors.  An  analysis  of  shear 
observations  in  an  isopycnal  or  isothermal  (constant  density 
or  temperature)  coordinate  frame  should  significantly  in¬ 
crease  our  understanding  of  shear  variability.  An  analysis 
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of  this  type  could  be  performed  using  existing  data.  For 
example,  the  MILE  data  set  discussed  in  this  report  would  be 
appropriate.  Shear  variability  due  to  vertical  displacements 
by  internal  waves  and  Ekman  pumping  could  be  largely  filtered 
out,  leaving  a  record  whose  variability  is  due  to  purely 
temporal  effects.  This  technique  would  be  ideal  for  studying 
the  autocorrelation  properties  of  shear,  and  the  coherence 
between  shear  and  wind  stress. 
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